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ABSTRACT 

We report on 10 years of monitoring of the 8.7-s Anomalous X-ray Pulsar 4U 0142+61 using 
the Rossi X-Ray Timing Explorer (RXTE). This pulsar exhibited stable rotation from 2000 March 
until 2006 February: the RMS phase residual for a spin-down model which includes is, v, and v 
is 2.3%. We report a possible phase-coherent timing solution valid over a 10-yr span extending 
back to March 1996. A glitch may have occured between 1998 and 2000, but is not required by 
the existing timing data. The pulse profile has been evolving since 2000. In particular, the dip 
of emission between its two peaks got shallower between 2002 and 2006, as if the profile were 
evolving back to its pre-2000 morphology, following an earlier event, which possibly also included 
the glitch suggested by the timing data. These profile variations are seen in the 2—4 keV band 
but not in 6—8 keV. We also detect a slow increase in the pulsed flux between 2002 May and 
2004 December, such that it has risen by 36±3% over 2.6 years in the 2—10 keV band. The 
pulsed flux variability and the narrow-band pulse profile changes present interesting challenges 
to aspects of the magnetar model. 

Subject headings: pulsars: individual(4U 0142+61) — stars: neutron — X-rays: stars — 



1. Introduction 

The existence of magnetars - young, isolated 
neutron stars powered by the decay of an ultrahigh 
magnetic field - is now well s upported by many in- 
depen dent lines of evidence (jWoods fc Thompson 
2006). This comes from the study of soft-gamma 
repeaters (SGRs) and anomalous X-ray pulsars 
(AXPs), both of which classically exhibit X- 
ray pulsations having luminosity in the range 
10 34-36 erg s _1 , periods P ranging from 5-12 s, 
Ps of 10 -13 — 10 -11 , and surface dipolar magnetic 
fields B in the range 0.6—7 x 10 14 G, assuming 
vacuum dipole magnetic braking. AXPs and 
SGRs, in the magnetar model, are ultimately 
powered by the internally decaying magnetic 
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field. In the magnetar model, the pulsed X-rays 
are suggested to be the result of a combination 
of surface thermal emission and a non-thermal 
high-energy component from resonant scattering 
of thermal photons off magnetospheric currents 
([Thompson. Lvutikov. fc Kul^rn J Eo0 2) . Magne- 
tar bursting, the hallmark of SGRs and also seen 
in AXPs, is believed to be a result of crustal yield 
and subsequent magnetospheric disturbances ulti- 
mately caused by stresses on the cr ust by the de- 



cayin g internal field (see for example I Woods et al 
200411 . 



Recently, thanks in large part to long-term 
monitoring campaigns, it has become clear that 
AXPs exhibit a variety of types of aperiodic X- 
ray variability that can in principle be useful for 
testing aspects of the magnetar model. This vari- 
ability can be categorized into four types, some 
of which are seen contemporaneously with each 
other: very short-duration SGR-like bursts, sud- 
den outbursts and transient brightenings with de- 
cays lasting months or longer, slow-rise long-term 
flux variations also with slow decays, and pulse 
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profile changes. 

Classic examples of SGR-like bursts and an out- 
burst were seen in 2002 for AXP IE 2259+586, 
which exhibited a sudden order-of-magnitude in- 
crease in the pulsed and tota l flux, followed b 



a one-year-long f lux decay (jKaspi et al.l 1200 



Woods et al.l 120041 ). The outburst was accompa- 



nied by over 80 short SGR-like bursts, a rotational 
glitch with interesting recovery on a time-scale of 2 
weeks, short-lived spectral changes, and dramatic 
broad-band pulse morphology changes which in- 
cluded the two profile peaks swapping heights and 
which lasted 2-3 weeks. The event was consistent 
with the picture of sudden crustal yield influenc- 
ing both the interior and the exterior of the AXP, 
in analogy with large SGR bursts. 

Several observations of this same source and 
others suggest similar outbursts in AXPs that 
went undetected. GI NGA observations of 
IE 2259+586 reported by|l wasawa. Kovama. fc Halpernl 
(|1992r ) also showed a factor-of-two pulsed flux 
change and pulse profile variations, both of which, 
in hindsight, could be explained by a contempora- 
neous but short-lived outburst that went unseen. 
In AXP 1RXS J170849.0-4000910, two rotational 



glitches were discovered flKaspi. Lackey, fc Chakr abartv 
20001: iKaspi fc Gavriill2003l:lDairOsso et al.ll2003f ). 



DaU'Osso et alJ ( 20031 ) reported possible small 



pulse morphology changes associated with these 
glitches. Whether these glitches were accompa- 
nied by bursting that went unobserved is unknown 
but plausible. The transient AXP XTE J1810-197 
underwent a dramatic sudden b rightening by 



nearl y two orders of magnitude ([Ibrahim et al 



2004) fol lowed by a total flux d ecay that lasted 
years (lHalpern fc GotthelJ l2005h : this may well 
have been an outburst similar to, though larger 
than that in IE 2259+586, but for which the brief 
main event went observed. Similarly, the transient 
candidate AXP AX J1845-0258 underwent a fac- 
tor > 100 decay in flux after an initial brighten 



ing that lead to its discovery ( Vasis ht et al 



ghten- 
■1 120001: 



Tarn et al.l 12009 ). This too could have been the 



result of an unseen outburst. 

AXP outbursts appear to be fundamentally dif- 
ferent from the slow-rise, long-ter m flux variations 
seen i n AXP IE 1048.1-5937. iGavriil fc Kaspi 
(|2004t ) discovered two long-lived, slow-rise X-ray 
pulsed flux flares from this source. The first flare 
had peak pulsed flux a factor of ~2 greater than 



the quiescent level, and lasted ^100 days. The 
second, larger flare had peak a factor of > 4 
higher than in quiescence, and lasted over one 
year. The flares were accompanied by an increase 
in the phase-averaged flux of the source and a de- 
crease in pulsed fraction, although the time scale 
and full dynamic range of the t hese changes have 
not b een clearly estab lished (jMereghetti et al 



2004 iTiengo et al.ll2005l ). No simultaneous pulse 



morphology changes were detected, and though 
the source did exhibit some SGR-like bursts 
( Gavriil. Kaspi. fc Woods! l2004t) . they were not 
obviously correlated with pulsed flux. Large (fac- 
tor of 10) torque changes were seen especially 
during the large flare, but the correlation between 
torque and pulsed flux, at least on time scales 
smaller than the flare itself, was marginal. Over- 
all, the slow-rise flares seen in IE 1048.1—5937 are 
not thought to result from crustal cracking as in 
outbursts. They can, however, be explained by a 
spontaneous increase in the magnetic field twist in 
the magnetosphere. However what mi ght trigger 
such events is unclear. Ne v erthel ess, iRea et al 



(|2005l ) and ICampana et al 
observations of 1RXS J170849.0 



(2007|) found, using 
400910 which 



has also shown phase-averaged flux variability, 
that one important prediction of the twisted mag- 
netosphere model appears to hold, namely a cor- 
relation between spectral hardness and flux. 

4U 0142+61 is an 8.7-s AXP. It has 
P = 0.2 x 10 -11 , implying a surface dipole mag- 
netic field of 1.3 x 10 14 G. From cont inuous RXTE 
monitoring, IGavriil fc Kaspil (|2002l ) showed that 



4U 0142+61 rotates with h i gh sta bility. However, 
Morii. Kawai. fc Shibazakl (|2005h reported a pos- 



sible timing glitch in 1999 on the basis of an ASCA 
observation in which the value of the frequency 
was margina lly discrepant w i th the f requency as 
reported by IGavriil fc Kaspil (|2002l ). iMorii et al. 
also reported simultaneous pulse morphology 
changes. As of 2006 April, in the published flux 
history of this source, there had been no reports of 
any X-ray activity like the flares of IE 1048-5937 
and the outburst of IE 2259+586. However, very 
recently, in 20 06 April and June and again in 
2007 Februar y dKaspi et alJliooj bib et alJliool 
Gavriil et al.l l2007tl SGR-like bursts were seen 



from 4U 0142+61, along with a sudden pulse pro- 
file change and a timing anomaly. In this paper 
we refer to the history prior to April 2006. 
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Here we report on continued RXTE monitor- 
ing observations in which we find a possibly new 
type of AXP variability, namely a slow, long-term 
increase in the pulsed flux accompanied by slow 
pulse profile evolution. We also report on f years 
of timing and investigate the possibility of a glitch 
having occured in f 998 or f 999, during an observ- 
ing gap, which may have precipitated the evolu- 
tion we are witnessing today. Our observations 
are described in Section [2 Our timing, pulsed 
morphology, and pulse flux analysis are presented, 
respectively, in Sections 13. 1[ 13. 2\ and 13.31 We 
present a combined pulse shape and pulsed flux 
analysis in Section 13.41 In Section [4j we discuss 
the possible origins of this behavior and the impli- 
cations for the magnetar model. 

2. Observations 

The results presented here were obtained us- 
ing the Proportional Counter Array (PCA) on 
board RXTE. The PCA consists of an array of 
five collimated xenon/methane multi- anode pro- 
portional counter units (PCUs) operating in the 
2—60 keV range, with a total effective area of ap- 
proximat ely 6500 cm 2 and a field of view of ~1° 



FWHM (jjahoda et al.lll996l ). Our 136 observa- 
tions are of various lengths (see Table [T]). Most 
were obtained over a period of several years as part 
of a long-term monitoring program, but some are 
isolated observations (see Fig. [T]) . 

Note that there is a 2-year gap in the observa- 
tions on which we are reporting: no RXTE obser- 
vations were made from 03 /21 /1998 (M JD 50893.083) 
to 03/07/2000 (MJD 51610.617). The gap exists 
because 4U 0142+61 was only added to our regular 
AXP monitoring program at the start of RXTE 
Cycle 5. Prior to the gap, our observations consist 
of a) 4 very closely spaced RXTE Cycle 1 observa- 
tions, b) 14 short Cycle 2 observations spanning 
a period of a year, and c) a single Cycle 3 ob- 
servation. No observations were made in Cycle 4 
(see Table H Fig. [J). 

For the monitoring, we used the GoodXenon- 
withPropane data mode except during Cycle 10 
when we used the GoodXenon mode. Both data 
modes record photon arrival times with l-/xs res- 
olution and bin energies into one of 256 channels. 
To maximize the signal-to-noise ratio, we analysed 
only those events from the top xenon layer of each 



PCU. 

3. Analysis and Results 
3.1. Phase-coherent Timing 

Photon arrival times at each epoch were ad- 
justed to the solar system bary c enter using the po- 
sition obtained by IPatel et all (|2003l ) from Chan- 
dra data. They were then binned with 31.25-ms 
time resolution. In the timing analysis presented 
below, we included only the events in the energy 
range 2—10 keV (unless otherwise specified) to 
maximize the signal-to-noise ratio of the pulse. 

Each barycentric binned time series was epoch- 
folded using an ephemeris determined iteratively 
by maintaining phase coherence; see below. Re- 
sulting pulse profiles, with 64 phase bins, were 
cross-correlated in the Fourier domain with a high 
signal-to-noise template created by adding phase- 
aligned profiles from all observations. The cross- 
correlation returned an average pulse time of ar- 
rival (TOA) for each observation corresponding to 
a fixed pulse phase. The pulse phase <f> at any time 
t can be expressed as a Taylor expansion, 



<P{t) = Mto) + Mt - to) + ^Mt ~ h) 



+-xv'o{t - t f + 
6 



(1) 



where v = 1/P is the pulse frequency, v = dv/dt, 
etc., and subscript "0" denotes a parameter eval- 
uated at the reference epoch t = t a . The TO As 
were fitted to the above polynomial using the pul- 
sar timing software package TEMPC0. 

We report an unambiguous phase-coherent tim- 
ing solution that spans the post-gap (i.e. after 
March 2000, MJD 51610) 6-yr period up until 
February 2006 (MJD 53787) including all data in 
RXTE Cycles 5—10. The parameters of our best- 
fit spin-down model which includes v, i>, and v are 
presented in Table [2 The corresponding phase 
residuals are shown in Figure [2 Note the unmod- 
elled features in the residuals; these may be caused 
by a noise process similar to th at commonly seen 
in ra dio pulsar timing (e.g. iLivingstone et al. 
120051 ). 

The best-fit post-gap ephemeris does not, how- 
ever, fit the pre-gap TOAs well. Figure [3] shows a 
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clear systematic deviation in the pre-gap residuals 
obtained after subtracting the post-gap ephemeris. 
The best-fit frequency obtained from the post- 
gap model at the reference epoch is larger than 
the frequency obtained from the best-fit model 
of the pre-gap TOAs at the same epoch (see Ta- 
ble [2|. This, in principle, could indicate that 
a glitch occured at some time during the gap. 
At MJD 51250, midway between the pre-gap 
and the post-gap ephemerides, the fractional 
change in frequency due to the possible glitch 
is Avjv = (7.11±0.15)xlCT 7 . However, by us- 
ing six frequency derivatives, we found a possible 
ephemeris that fits the entire Cycle 1 to 10 range 
(MJDs 50170 to 53787, see Table H}. The RMS 
phase residual for that ephemeris is 0.019 (see 
Fig. [4|. Note that when finding an ephemeris that 
spans several years, it is not uncommon to require 
a large number of frequency derivatives in order to 
reduce the RMS phase residuals to a number on 
the order of 5%. This is typical especially in young 
puls ars and is generally attri buted to timing noise 



puls ars ana is g enerally attrib 
(e.g. iLivingstone et al.ll2005h . 



The possibility of a glitch in 4U 142+61 dur- 
i ng ou r gap was in fact examined by M orii et al 
(|2005h . The authors showed that the frequency 



obtained from a 1998 August ASCA observation 
of 4U 0142+61 (MJD 51046.699; 154 days after 
the start of our gap), differs from the frequency 
predicted at the epoch of the observation by the 
e phemerides previously r e ported for 4 U 0142+61 



iicmcricies previously r e ported lor 4 U U14Z+0 
Gavriil fc Kaspl (|2002h . iMorii et all reported ; 



frequency / = 0.1150972(6) Hz at MJD 51046.69875 
for the ASCA observation. Our overall ephemeris 
(see Table EJ) predicts / = 0.115098404(3) Hz at 
the same epoch. Their measurement is within 
2a of our prediction, indicating a ^5% possibility 
that the two values are the same. Therefore, their 
measured / can be explained by a gradual change 
of the spin-down rate without invoking a glitch. 

However, the existence of our overall ephemeris 
cannot rule out the possibility of the glitch: in 
some rotational glitches, the frequency evolution, 
given some relaxation time after the glitch epoch, 
returns to what it was prior t o the glitch (see 
for ex ample the glitch reported in lKaspi fc Gavriij 
2003) . If a glitch of this kind had happened inside 
the two-year gap, and if the length of the gap was 
much greater than the relaxation time, the only 
long-term effect of the glitch that could still be ob- 



servable with a timing analysis would be a random 
phase jump in the post-relaxation TOAs relative 
to the pre-glitch TOAs. To investigate this pos- 
sibility, we added an arbitrary but constant time 
jump to all the post-gap TOAs. We were still able 
to find a new ephemeris that connected the TOAs 
through the two-year gap. This indicates that our 
overall ephemeris is not unique. Hence, we cannot 
rule out the possibility of a random phase jump 
between Cycles 3 and 5, and therefore, a glitch 
cannot be ruled out. 

It is important to note that our method for ob- 
taining TOAs, (cross-correlating the folded pro- 
files of given observations with a high signal-to- 
noise template obtained from all the observations 
combined), assumes a constant pulse profile. In 
the next Section, we show that the pulse profile is 
actually varying during our monitoring program. 
However, we performed simulations which showed 
that these changes do not result in timing offsets 
significantly larger than the reported TOA uncer- 
tainties. Hence the profile variations do not affect 
the above analysis. 

3.2. Pulse Profile Changes 

3.2.1. Qualitative Observations 

We performed a pulse profile analysis using 
FTOOLS version 5.3. 10. We used the following 
steps: for each observation, we ran the FTOOL 
make_se to combine the GoodXenon files. We 
then used the FTOOL f asebin to make a phase- 
resolved spectrum of the entire observation with 
64 phase bins across the profile. When we ran 
f asebin, we selected layer 1 of the detector, dis- 
regarded the propane photons, and included the 
photons from PCUS 1, 2, 3, and 4. We disre- 
garded photons from PCU because of the loss 
of its propane layer in 2000 (jjahoda et al. 20061 ) 
and because it gave different results from the other 
PCUs. f asebin also took care of barycentering 
the data. For each observation, we then used 
seextrct to make a phase-averaged spectrum for 
the same set of detector layers and PCUs. The 
phase-averaged spectrum was then used by the 
perl script pcarsp to make a response matrix. 

We loaded the phase-resolved spectra and the 
response matrices into the X-ray Spectral Fit- 



3 http: / /heasarc. gsfc.nasa.gov/ftools 
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ting Package (XSPEC0) and selected photons 
belonging to three energy bands: 2—10, 2—4, 
and 6—8 keV. Using XSPEC, we extracted an 
ASCII count rate pulse profile for each of the 
energy bands. The profiles included XSPEC- 
obtained la error bars on each of the phase bins 
in the profiles. To obtain a pulse profile in units of 
count rate per PCU, we divided the overall profile 
by a PCU coverage factor that took into account 
the amount of time each PCU was on. 

We then aligned the 64-bin profiles with a 
high signal-to-noise template using a similar cross- 
correlation procedure to the one described in the 
timing analysis. Then, for each RXTE Cycle, we 
summed the aligned profiles, extracted the DC 
component from the summed profile, and scaled 
the resulting profile so that the value of the high- 
est bin is unity and the lowest point is zero. 

The average profiles in all three bands are pre- 
sented in Figure[5]for comparison. In a given band, 
the different profile qualities are due to different 
net exposure times. 

It is important to note that the two narrow 
energy bands that we are using contain photons 
belonging to different spectral components: from 
the spectrum of 4U 0142+61 (see, for example, 
White et al.lll996h . under the assumption that the 



spectrum is well described by a blackbody plus 
power-law tail, we know that the higher energy 
band (6—8 keV) contains negligible blackbody 
emission, while the lower energy band (2—4 keV) 
contains comparable amounts of blackbody and 
power-law emission. 

Qualitatively, the evolution of the pulse profiles 
in the first two bands of Figure [5] is clear to the 
eye. In Cycles 1 and 2, the smaller peak, obvious 
in later Cycles, is not very well defined. After the 
two-year gap, in Cycle 5, the 'dip' between the 
peaks is much more pronounced. The emission in 
the dip starts to most noticeably rise in Cycles 
subsequent to Cycle 7. In the 6—8 keV band, the 
smaller peak appears to have lower amplitude in 
the normalized profiles than in 2—4 keV, indicat- 
ing that it has a softer spectrum relative to the 
larger peak. 

Another qualitative observation is that the ra- 
tio of the heights of the two peaks in the 2—10 keV 
and 2—4 keV bands appears to be closest to unity 



4 http://xspec.gsfc. nasa.gov Version: 11.3.1 



in the first Cycle after the gap. The ratio starts 
to decrease in the Cycles subsequent to Cycle 5. 
Note that from this Figure alone we can compare 
the sizes of the two peaks, but we cannot track the 
evolution of the heights of each peak separately. In 
order to do that, we need to scale the pulse pro- 
file of each Cycle by the average pulsed flux. This 
analysis is presented in Section [3.41 

3.2.2. Fourier Analysis 

To quantify the changes in the pulse profile, wc 
computed the first six Fourier amplitudes of the 
average profiles of each Cycle in each energy band. 
Harmonic numbers larger than 6 were always con- 
sistent with zero. The results are shown in Fig- 
ures El and[H In each of the three Figures, the 
plots on the right show the power in each harmonic 
divided by the total power in all harmonics (not 
including the DC term). In the plots on the left, 
the observed pulse profiles are in the background 
(hollow squares without error bars). The super- 
imposed smooth curve in the foreground is made 
from the first six calculated Fourier components. 
The ratios of the first three Fourier harmonics rel- 
ative to the fundamental are presented in Figure[9] 
Note that the ratios of the Fourier harmonics are 
not presented for Cycle 3 in Figures [TJ and [8] 
due to the very low signal-to-noise ratio. 

Significant variations were seen in the pulse pro- 
file of 4U 0142+61 in the 2-10 keV band and in 
2-4 keV; on the left side of both Figures [6] and [3 
the most striking variable feature is the difference 
in the relative heights between the top of either 
peak and the bottom of the dip. On the right side, 
considering the first two Fourier components, the 
ratio of the second to first amplitudes (A2/A1) is 
significantly bigger than unity only in Cycle 5. It 
then falls steadily until, in Cycle 10, it reaches the 
same ratio as in Cycle 1. The evolution of the 
A2/A1 ratio in 2—4 keV is shown in Figure [5] In 
the pre-gap Cycles, harmonics of order higher than 
2 are only marginally significant. In the post-gap 
Cycles, harmonics 3 and 4 are most significant in 
Cycles 5 and 6, coinciding with the Cycles where 
the dip in the time-domain curve is sharpest. The 
evolution of the A3/A1 and the A4/A1 ratios is 
also shown in Figure [H Note the obvious rise 
in the harmonic ratios just post-gap in the soft- 
band, with apparent subsequent evolution to pre- 
gap values. 
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In the 6—8 keV band (see Fig. [8]), statements 
about the behavior of the Fourier components are 
harder to make because of the poor signal-to-noise 
ratio; nevertheless some trends are clear. Unlike 
in the lower energy band, the A 2 /Ai ratio does 
not appear to systematically increase or decrease. 
Also unlike in the lower band, harmonics 3 and 4 
do not appear to vary systematically (see Fig. [9]). 
Thus, in the band where all the emission is from 
the power-law component of the spectrum, the 
variations in the shape of the pulse profile, if any, 
are much less significant than the variation in the 
2—4 keV band which contains photons belonging 
to both components of the spectrum. 

3.3. Pulsed Flux Time Series 

To obtain a pulsed flux time series for 4U 0142+61, 
we did the following. First, for each PCU in each 
observation, we used a procedure similar to that 
described in Section 13.2.11 to make a phase re- 
solved spectrum (with 16 phase bins across the 
profile) and a response matrix. We then used the 
FTOOL f modtab to correct the exposure value in 
the phase resolved spectrum of each PCU in or- 
der to take into account the amount of time that 
each PCU was on. Then, for each observation, 
we added the spectra obtained from PCUs 1 to 4 
using the FTOOL f badd and added the responses 
using the FTOOL addrmf . We used fbadd and 
addrmf again to add the spectra and responses 
of all observations in a given RXTE Cycle. For 
each RXTE Cycle, we loaded the phase resolved 
spectra into XSPEC, and selected photons in the 
2—10 keV range. Using XSPEC, we extracted 
an ASCII count rate pulse profile for each RXTE 
Cycle. The profiles included XSPEC-determined 
la error bars on each of the phase bins. We then 
smoothed each of the profiles by eliminating the 
Fourier components corresponding to harmonic 
numbers larger than five. The pulsed flux for each 
of the smoothed profiles was calculated using the 
following discrete area formula: 



F 



N 



i=l 



{jPi Pmin) • 



(2) 



where i refers to the phase bin, iV=16 is the total 
number of phase bins, Pi is the count rate in the 
i^h phase bin of the smoothed pulse profile, and 
Pmin is the value of the minimum of the continuous 



smooth function that is made from the first five 
Fourier components of the original profile. 

The resulting pulsed flux history in counts/s/PCU 
is shown in the top panel of Figure [TO] Each 
point represents one RXTE Cycle. The pulsed 
flux has increased by 36±3% between Cycles 7 
and 9. A quick rebinning of the observations 
shows that the increase period lasted ~ 2.6 yr 
(between MJDs 52400 and 53350). We verified 
that the same trend is detected in PCUs 1—4 indi- 
vidually, and that there are no comparable trends 
in the long-term light curves of the other 4 AXPs 
observed as part of this monitoring program. 

We repeated the above procedure of finding 
the flux for narrower energy bands. There are 
hints that the long-term increase is present in the 
2—4 keV band and not in 6—8 keV but our statis- 
tics do not let us confirm this. If the pulsed flux 
increase is restricted to < 6 keV, this could in- 
dicate that the spectrum of the pulsed emission 
is getting softer. Motivated by this possibility, 
we performed a detailed spectral analysis of four 
available archival XMM observations. We found 
that the spectrum is indeed getting softer (Gon- 
zalez et al. in preparation). 

Note that the method that we used to calcu- 
late the pulsed flux, which consisted of calculating 
the pulsed area under the profile, is more sensi- 
tive to noise than are measurements of the RMS 



pulsed flux like those used in I Woods et al.1 (|2004[ ). 
Therefore, to reduce the effects of noise, it was 
necessary to combine the data from entire Cycles 
in order to obtain each of the pulsed area points 
reported in the top panel of Figure [HJ hence the 
large horizontal error bars. We report measure- 
ments of the pulsed area instead of RMS pulsed 
flux because, while it is true that measurements 
of the latter are less sensitive to noise, changes in 
it can be caused by changes in the real pulsed area 
and also by variations in the pulse profile. As a 
double check, we calculated the RMS pulsed flux 
for each observation, found the average of the re- 
sulting fluxes in each RXTE Cycle, and multiplied 
each average by a conversion factor dependent on 
the pulse shape, in order to obtain a pulsed area. 
When we plotted these pulsed areas versus time, 
we obtained a similar trend to that see in the top 
panel of Figure [TO] For further discussion of the 
various methods used to estimate the pulsed flux 
see Archibald et al. (2007, in preparation). 
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In order to verify that the trend seen in the 
top panel of Figure [10] is not an artifact of the re- 
sponse of the detector, and in order to verify that 
the relative levels in the pre-gap and the post- 
gap flux values are not skewed due to an evolu- 
tion of the response of the detector, we calculated 
the pulsed flux in erg/s/cm 2 using an additional 
method that took the evolution of the response 
into account. For each RXTE Cycle, we obtained 
one phase resolved spectrum (with 16 phase bins 
across the profile) and one response matrix using 
the method described above. This time we in- 
cluded PCU 0. Then, we defined the background 
phase bin to be that where p m in h es - We then used 
the FTOOL cmppha to combine the spectra of all 
phase bins into a single phase-averaged spectrum. 
Then, using XSPEC, we subtracted the spectrum 
of the background bin from the combined spec- 
trum of the remaining phase bins, scaling the ex- 
posure appropriately. When subtracting the spec- 
trum of the background bin from the spectrum of 
the remaining bins, we assumed that the result- 
ing spectrum is, to a good approximation, that of 
the pulsed component of the emission. For each 
RXTE Cycle, we fit the spectrum of the pulsed 
emission with a model that consisted of a black 
body plus a power law. When doing the fitting 
we froze the temperature of the black body to the 
value kT = 0.44 keV, the power-law photon in- 
dex to 7 = 3.94 and the column density of neutral 
hydrogen to Nh — 0.99 xlO 22 cm -2 . These pa- 
rameters were obtained from a linked spectral fit 
of the four archival XMM observations (see Gon- 
zalez et al. in preparation). We let the normal- 
izations of the two components of the spectrum 
vary freely. After the fit was done, we extracted 
the pulsed flux numbers from the fitted spectrum 
of the pulsed emission of each RXTE Cycle. We 
multiplied the la error bars returned by XSPEC 
by the square root of the reduced chi squared of 
the spectral fit. This multiplication is equivalent 
to assuming that the model being fit is the right 
one to use and that any bad-fit results from an ini- 
tial underestimation of the error bars of the spec- 
trum that is being fitted. The resulting pulsed flux 
values in erg/s/cm 2 for each of the RXTE Cycles 
are presented in the bottom panel of Figure [TO] 
The trend observed in the pulsed flux is similar 
to that in the top panel. The observed increase 
between Cycles 7 and 9 is 29±8%, consistent with 



that found in the first analysis. However, after 
having taken the response of the detector into ac- 
count, the pulsed flux in the pre-gap observations 
appears to be consistent with the pulsed flux in 
the first post-gap Cycle. 

From the bottom panel of Figure [TO] the in- 
crease in the pulsed flux is ~ 2.1 x 10~ 13 erg s" 1 cm" 
in the 2—1 keV band. Assuming a distance 
of 2. 5 kpc ( Hulleman. van Kerkwiik. fc Kulkarnil 
20041 ), the total luminosity increase in the ~ 2.6- 
yr period during which the increase happened 
is ~ 1.1 xlO 33 erg s _1 . This increase is of the 
same order of magnitude as the average en- 
ergy release rate in th e first IE 1048—5937 flare 
(|Gavriil fc Kaspil 120041) . It is also an order of 
magnitude smaller than the average energy re- 
lease rate in the first day following the outburst 
in IE 2259+586 | Woods et al.ll2004f l. The amount 
of energy released in the same 2.6-year period due 
to the increase in the pulsed flux is ~ 9xl0 40 erg 
in the 2—10 keV band. This is comparable to 
the e nergy released in the s econd IE 1048—5937 
flare ( Gavriil fc Kaspil 2004h . and to the energy 
released in the i n the year following the outburst 



in IE 2259+586 ( Woods et alJ l20041 



3.4. 



Combined Pulse Morphology and 
Pulsed Flux Analysis 



In Section [3. 2. 21 we calculated the Fourier com- 
ponents of the average pulse profiles. This gave us 
the relative amplitude of the pulse profile harmon- 
ics in each RXTE Cycle. In Section l3~3] we calcu- 
lated the pulsed flux for every observation. Here, 
we compute a weighted average of the pulsed flux 
for each Cycle using the flux points calculated in 
Section 13.31 We then reconstruct the profiles for 
each of the Cycles from the first six Fourier com- 
ponents (not including the DC), scale them by the 
average RMS pulsed flux for that Cycle, and add 
the necessary offset for the lowest point on each 
curve to be zero. This means that the resulting 
scaled profiles return the correct pulsed flux. The 
advantage of this analysis is that we can now trace 
the evolution of each of the peaks independently. 
The post-gap scaled profiles in 2—10 keV and in 
2—4 keV are presented in the top panels of Fig- 
ures [TT] and [TOl respectively. We did not include 
a similar Figure for 6—8 keV because of the poor 
signal-to-noise ratio in that band. The absolute 
heights of the peaks in the post-gap Cycles, as well 
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as the absolute height of the dip in between, are 
plotted in the bottom panels of Figures fTTI and fT2l 
The error bars take into account both the errors 
on the Fourier components and the errors on the 
pulsed flux. 

In both Figures, there is a hint of increase in 
the height of the big peak between Cycles 7 and 9. 
The dip between the peaks appears to be getting 
shallower more rapidly between Cycles 7 and 9. 
The difference in the height of the dip over these 
3 years is more significant than the difference in 
the height of either peaks. This indicates that the 
biggest contribution to the change in the pulsed 
flux comes from an increase in the emission in the 
dip, which, in principle, could be caused by the 
widening of either peak around the dip. 

4. Discussion 

4.1. Possible Event in the Gap? 

Could a short-time scale energetic event (such 
as an outburst like that seen in 2002 for IE 2259+586) 
have occured sometime within the two-year gap 
and triggered the pulsed flux and pulse profile 
changes that we are observing? As discussed 
above, the possibi lity of a glitch dur ing the gap 
was examined by M orii et al.l (|2005h . Here, an 
examination of our timing, flux, and pulse profile 
analyses can provide further clues to help answer 
this question. 

From our timing analysis (Section [3TT]), there is 
some evidence for a glitch having occured some- 
time during our gap. Hence, if an outburst did oc- 
cur, it might have been accompanied by a glitch, as 
was the case for the 2002 outburst of IE 2259+586 
(jKaspi et al. I l2003h . If there was a glitch in 
4U 0142+61, the unrecovered fractional change 
in frequency would have been (7.11+0.15) x 10 -7 , 
a factor of 6 smaller than the maximum fractional 
frequency change of (4.24+0.11) xlO -6 observed 
in IE 2259+586. The fact that the pulsed flux in 
the first post-gap observations is consistent with 
that in the last pre-gap observations could be con- 
sistent with an outburst in between if the initial 
flux increase during an outburst had time to die 
down (see Fig. [TTJ]). If we assume that the re- 
turn of the pulse profile to its pre-gap shape is 
a recovery following an outburst, this would im- 
ply a much longer timescale for the pulse profile 
relaxation phase than for the pulsed flux relax- 



ation, the opposite to what was seen fo l lowin g 



the IE 2259+586 outburst (| Woods et al J 12004 - 



Alternatively, the post-outburst pulse profile re- 
laxation could have been completed during the 
gap, and the slow return of the post-gap profile 
to its pre-gap morphology could be attributed 
to a different phenomenon. This is further dis- 
cussed in Section 14.31 In either case, if there was 
an event during the gap, why the pulsed flux is 
presently rising is unclear. If the event associated 
with the putative glitch released energy deep in 
the neutron-star crust, the n the increase could be 
due to its slow rele ase (eg. Eichler fc Cheng 1989t 
Hirano et al ] ll997h . Given the size of the observed 



flux increase and its timescale for release, the ini- 
tial energy deposition would have had t o have 
been large, ~10 45 erg (jHirano et al.lll997l ). This 
is comparable to the observed ener gy release in 
giant SGR flares (|Hurlev et al.ll2005l ). 

If we assume that the pulse shape prior to the 
gap is the "relaxed" pulse shape, the evolution of 
the harmonic ratios shown in Figure [9] supports 
the possibility of relaxation of the profile following 
an event in the gap. To shed light on the events in 
the gap, we can compare our RXTE profiles wit h 
those observed with ASCA bv lMorii etaH (|2005l ). 
In their Figure 4, pulse profiles in 0.5—10 keV for 
a) September 1994, b) August 1998, and c) com- 
bined July and August 1999 profiles are presented. 
In the 1994 and 1998 observations, the profile con- 
sisted of two peaks, with the trailing peak being 
the smallest, with the dip between the peaks be- 
ing higher than the lowest bin in the profiles. The 
shape of the 1998 profile is in agreement with the 
RXTE pre-gap average pulse profiles for Cycle 2 
(see Fig. [B]). In the 1999 ASCA profile, the ampli- 
tude of the trailing peak was higher than that of 
the leading peak. In addition, the difference be- 
tween the height of the dip and the lowest point 
in the profile decreased. Interestingly, the changes 
in the ASCA profiles appeared more significant at 
the lower end of the energy band, as we observe 
in our RXTE data. In 2000, the first RXTE Cy- 
cle after the gap has a profile in which the trailing 
peak is once again smaller than the leading peak. 
The dip between the peaks, however, is still more 
pronounced than in the pre-gap observations. 

Overall, the timing data and the pulse profile 
data are consistent with some sort of event, possi- 
bly a glitch with accompanying sudden pulse pro- 
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file change, having occured between 1998 August 
and 1999 July, possibly with the latter's long-term 
relaxation still ongoing as of early 2006. However, 
we suggest an alternate explanation for the latter 
point below. 

4.2. Brief Review of the Magnetar Model 

In the detaile d mag netar model proposed by 
Thompson et al.l ( 2002f ). the crust of a magne- 
tar is deformed by internal magnetic stresses, 
thereby twisting the footpoints of the exter- 
nal magnetic field, driving powerful currents in 
the magnetosphere and twisting the magneto- 
sphere relative to the standard dipolar geome- 
try. These magnetospheric currents resonantly 
cyclotron scatter seed surface thermal photons. 
The seed contribution to the thermal component 
of the spectrum is thought to arise from heat 
resulting from the a ctive decay of a high inter- 
nal magnetic field (iThompson &; Duncan 1996; 
Thompson et alJl2002 )T The magnetospheric scat- 
tering is responsible for the non-thermal compo- 
nent of AXP spectra. Additionally, the surface 
is back-heated by the currents, resulting in ad- 
ditional thermal emission. Indeed, the persistent 
emission in AXPs generally has a spectrum that 
is well described by a two-component model, con- 
sisting of a blackbody plus a hard power-law tail, 
as expected in this model. 

Changes in pulsed and/or total X-ray lumi- 
nosity, spectral hardness, and torque are pre- 
dicted to have a comm on physical origin in the 
Thompson et~aT1 (2002) model and some correla- 
tions are expected. Changes in twist angle of 
the magnetic field, cause, or may be caused by, 
changes in the magnetospheric current distribu- 
tion (due either to sudden crustal deformation like 
in AXP outbursts and SGR giant flares or due to 
slower crustal defo rmations as may be tak ing place 
in the AXP flares; iGavriil fc KaspHl2004l) . Larger 
twists generally correspond to harder persistent 
X-ray spectra and higher magnitudes of the spin- 
down rates, as is observed when comparing the 
harder SGR spectra to t hose of the softer AXPs 
( Marsden fc White] 2001 ). A similar trend might 
be expected for a single magnetar exhibiting lumi- 
nosity variations: a higher luminosity should cor- 
respond to a larger twist, hence harder spectrum, 
as has been reported for 1RXS J170849 .0-400910 
(|Rea et al.l 120051: ICampana et al.ll2007t) . A higher 



luminosity should also in general correspond to a 
larger magnitude of the spin-down rate. However, 
decoupling between the torque and the luminos- 
ity can occur because the torque is most sensitive 
to the current flowing on a relatively narrow bun- 
dle of field lines that are anchored close to the 
magnetic pole. For a single source, whether an X- 
ray luminosity change will be accompanied by a 
torque change depends on where in relation to the 
magnetic pole the source of the enhanced X-rays 
sits. 

The Thompson et al. (2002) model can also ex- 
plain properties of the pulse profiles of magnetars. 
According to the model, several effects can affect 
the pulse shape, generate subpulses, and/or in- 
crease the energy dependence of the pulse profile. 
In addition, an increase in the twisting angle of the 
magnetic field increases multiple scattering and in- 
creases the optical depth to resonant scattering 
which can simplify the pulse shape. This is one 
of the two proposed explanations for the sudden 
simplification of the pulse pro file of SGR 1900+1 4 
after its dramatic giant-flare ([Woods et al.ll200lh . 
the other explanation being the sudden elimina- 
tion of the nonaxisymmet ric components of th e 
magnetospheric currents (jThompson et al.ll2002l ). 



4.3. 



Possible Physical Interpretations for 
4U 0142+61 



In this paper, we have shown that the pulsed 
flux of 4U 0142+61 has increased on a time scale 
of a few years, and we have found simultaneous 
slow pulse profile evolution in the 2—4 kcV band, 
which may to be recovering from some event that 
occured prior to 2000 but after 1997, possibly in 
the interval between 1998 August and 1999 July. 
We have als o foun d evidence, as first suggested by 
Morii et al.l (|2005l ). that there was a timing glitch 
in that same interval, although we cannot confirm 
its existence. Can the magnetar model explain 
these observations? 

The energy dependence of the pulse profile evo- 
lution is puzzling. As describ ed above, in the 



twiste d magnetosphere model of iThompson et al 
(|2002h . the non-thermal emission in an AXP is 



the result of magnetospheric scattering of surface 
thermal photons. If the surface emission angular 
pattern were changing, therefore, the non-thermal 
angular pattern should as well. That we do not 
observe comparable pulse profile changes in the 
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6-8 keV band for 4U 0142+61 is thus puzzling. 
One possibility is that the seed thermal emission is 
not changing appreciably, but the scattering cur- 
rents in the outer regions of the magnetosphere, 
where the cyclotron energy is lower, are chang- 
ing, while the inner currents are not. This could 
arise if there is evolution in the field configuration 
closer to the magnetic poles, with relatively little 
closer to the magnetic equator. Why variations 
in the field configuration should be geographically 
localized, however, is unclear. 

The increase in the pulsed flux over a similar 
time scale as that of the profile evolution is appar- 
ently accompanied by a softening of the spectrum 
(Gonzalez et al. in preparation). As discussed 
above, the putative 1998/1999 glitch may have de- 
posited a large amount of energy in the crust, with 
it only starting to be radiated away in 2002. The 
energy rele ased would have ha d to have been large, 
~10 45 erg ( Hirano et al.lll997 ) . Such a thermal en- 
ergy release could be influencing the pulse profile 
as well, although some change in profile would be 
expected in the hard band too, which is not ob- 
served. Moreover, such an increase has not been 
observed following the glitches in IE 2259+586 or 
1RXS 170849.0-400910, although this could be a 
result of smaller total energy releases. 

Alternatively, the increase in pulsed flux seen 
in 4U 0142+61 between 2002 and 2005 could be 
explained by the twisted magnetosphere model. In 
this framework, there are two possibilities: a) a 
slow increase in the twist of the magnetic field 
lines in the magnetosphere, or, b) a slow decrease 
in the twist angle. 

In the first possibility, the observed increase in 
the pulsed flux could be an extreme case of the 
IE 1048-5937 flares, i.e. a slow twisting of the 
magne tic field lines in t he su ggested interpreta- 
tion of iGavriil fc Kaspi ( 2004 ). This explanation 
is only valid if the total flux, which remains to 
be determined with an imaging telescope, is in- 
creasing as well. How this would be related to the 
putative 1998/1999 event is unclear; any flux en- 
hancement that occured then would have had to 
have largely decayed away by 2000. In any case, 
for an increase in t he twist angle, one exp ects a 
spectral hardening (jThompson et al. 2002). We 
do not observe this. Also for an increase in the 
twist angle, at least naively, the twisted magneto- 
sphere model predicts an increase in torque as the 



flux rises (jThompson et al.ll2002h . From Table [2j 
the post-gap ephemeris v is positive, meaning v 
is increasing, i.e. the magnitude of the pulsar's 
spin down rate is decreasing, the opposite of what 
is expected for a flux increase, unless the mag- 
netospheric currents causing the torque are flow- 
ing only in the small polar cap region. Finally, 
if the slow increase in pulsed flux were caused by 
a slow magnetospheric twisting, the decrease in 
the size of the Fourier components of order higher 
than unity in the low energy band could be in- 
terpreted as a simplification of the pulse profile 
due to an increase in the scattering in the mag- 
netosphere. However, the most extreme case of a 
pulse profile simplification, which was reported in 
SGR 1900+14, happened equally in all bands fol- 
lowing a dra matic increase in the scattering after 
a gia nt flare ( Woods et al.ll200ll iThompson et al. 
2002). The phenomenon that we have observed, 
by contrast, is restricted to the softer energies. 

In the previous paragraph, we mentioned that 
the increase in the pulsed flux and the gradual 
changes in the pulse profile may indicate stress 
build-up caused by an increase in the twist an- 
gle in the magnetosphere. The pulse profiles of 
Cycles 2 and 10 are very similar, indicating that 
the pulse profile of Cycle 2 may also be showing 
signs of the same kind of stress build-up. Un- 
der these assumptions, if an event occured in the 
gap following Cycle 2, it is reasonable to expect 
a similar event to follow after Cycle 10. Indeed, 
in 2006 April, less than 2 months after the end 
of Cycle 10, the pulsar appears to have entered 
an extended active phase: a single burst accom- 
panied by a pulse profile change was detected 
from the pulsar on April 6 (jKaspi. Dib. fc Gavrii] 
2006). A series of four bursts was later detected 



on June 25 (|Dib et alj [2006) and a larger burst 



was d etected on 2007 February 07 (jGavriil et al 



2007). A detailed paper on these events is cur- 



rently in preparation. If they are indeed due to a 
stress release following several years of slow mag- 
netospheric twist, then the causes of the softening 
of the spectrum and of the decrease in the magni- 
tude of the pulsar's spin-down are unclear. 

The second possibility in the framework of the 
twisted magnetosphere model, is that the observed 
increase in the pulsed flux is accompanying a 
slow decrease in the twist angle of the magnetic 
field lines in the magnetosphere. This explana- 
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tion is only valid if the total flux is decreasing. 
If the total flux were falling with the pulsed frac- 
tion rising, the observed spectral softening could 
be consistent with naive untwisting expectations 
as would the decreasing spin-down rate. Indeed 
an anti-correlation betwe en total flux and pu lsed 
fraction was observed by iTiengo et al.1 (|2005l ) for 
IE 1048.1—5937, supporting this possibility. How- 
ever, the problem of why the pulse profile is evolv- 
ing at low, but not high energies, remains. 

4.4. Other Wavelengths 

4U 0142+61 is truly a multi-wavelength AXP. 
It is known to pulsate in the optica l band 



(Kern & Martin 2002: Dhillon et al. 2005) and it 



has been detected in the near-IR ( Hullcman et al 



2004), in the mid-IR using SPITZER (Wang, Chakrabarty 



2006) , and in hard X -rays (|Kuiper et all 12006: 



den Hartog et al.| [2006) 



The origin of the emission at these other 
wavelengths remains unclear, although some 
models have been proposed. I n hard X-rays, 
Thompson fc Beloborodov ( 20051 ) argue that the 
emission is either due to bremsstrahlung photons 
emitted by a thin surface layer, or is due to syn- 
chrotron emission originating from the region in 
the magnetosphere where the electron cyclotron 
energy is in the keV range. The near-IR and op- 
tical emission is thought to be mag netospheric 
(|Eichler. Gedalin. fc Lvubarskvl 120021 ) . while the 
mid-IR emission is suggested to be d ue to a pas- 
sive fall-back disk (jWang et al. I l2006l ). 

Looking for correlations between the X-rays 
and the emission in other wavelengths may serve 
as tests of emission models, as correlations be- 
tween X-ray and near-IR fluxes have sometimes 
been observed. For example, a correlation in 
the decays of the X-ray and near-IR fluxes in 
IE 2259+ 586 was observed following the 2002 
outburst ( Tarn et alj 12004 ): a similar correla- 



tion w as reported for XTE J1810-197 (|Rea et al 
20041 ). However, in other instances 



fluxes were not correlated, as for IE 



the two 
1048-5937 



Wang et al. (|2006h observed mid-IR emission 
from 4U 0142+61 which they argue is associated 
with a passive fall-back disk irradiated by the cen- 
tral X-ray pulsar. If this is the case, then if 
4U 0142+61 's X-ray flux is increasing, one expects 
a corresponding increase in the disk emission. It 
is thus important to establish the behavior of the 
total flux, in addition to that of the pulsed flux 
reported on here. 

The possible presence of a disk suggests that if 
a sudden, impulsive outburst occured in the gap, 
the energy released must have been significantly 
smaller than the disk binding energy, after ac- 
counting for the disk thickness. For a central pul- 
sar mass of Mpsr ~ 1-4M Q , a uniform disk of mass 
M ~ 3M e , and inner and outer r adii R\ ~ 3R Q 
agdpgj){>^lOR0 ( Wang et"ai1l2006h the disk bind- 
ing energy is ~4xl0 42 erg. Since the X-ray lumi- 
nosity of the source (^10 35 erg/s) integrated over 
a period > 5000 yr is much larger than the bind- 
ing energy of the observed disk, one must assume 
that when the source is not undergoing an out- 
burst, the disk is in an equilibrium state where the 
rate of energy absorption is balanced by the rate 
of disk emission. Assuming this equilibrium can- 
not hold on the timescale of a sudden outburst, 
then IR observations of the disk provide an up- 
per limit of (4xl0 42 //) erg on the energy released 
in a possible outburst in the gap, where / is the 
fraction of the solid angle occupied by the thick- 
ness of the disk. For / = 0.01, this upper limit 
is three orders of magnitude larger than the total 
energy released during the flares of IE 1048—5937 
(lOavriil fc Kasnill2004h . It is also five orders of 
magnitude larger than the energy released dur- 
i ng the first day o f the IE 2259+586 outburst 
(| Woods et al.l I2004T ) . Thus, an event in the gap 
of either the magnitude of the flares or that of the 
outburst occuring could have affected the disk but 
seems unlikely to have disrupted it. For / = 0.01, 
the upper limit is also of the same magnitude as 
the energy released in either of the SGR 1900+14 



(iGayriil fc Kasoil I2004L iMereghetti et all I2004L 
and Durant. van Kerkwiik. fc Hullemanll2004h . In 
4U 0142+61, reported variations in the IR are 
not seen contemporaneously in our X-ray data, 
and are on time scales much shorter than that of 



or th e SGR 0526—66 giant flares (|Mazets et al 
1999). This suggests that there have been no 



the X-ray variation reported here (|Hulleman et al 
20041 iDurant fc van Kerkwijkll2006l )~ 



events of the magnitude of the giant SGR flares 
since the putative disk's formation. Given that 
we have witnessed 3 giant SGR flares each from a 
different source since 1979, and none from AXPs, 
this suggests that AXPs do not exhibit giant flares. 
This would suggest an interesting distinction be- 
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tween SGRs and AXPs: if the two are evolutionary 
linked, the AXP phase must come first, unless a 
debris disk can reform following a giant flare. 

The mid-IR emission has also been interpreted 
as an active fall-back disk in which the pulsar is 
accreting in a propeller mode. In this case, the 
dipole field strength of the pulsar is ty pical of con- 
venti onal radio pulsars, ie. ~ 10 12 G (jErtan et al 
2007), although the surface field strength is in 



the magnetar range, owing to higher order mul- 
tipoles. In this model, the X-rays arise from pro- 
peller accretion. Without detailed models of the 
spectra and of the pulse shapes expected in this 
model, we cannot interpret our observations in 
this framework. Nevertheless, one might expect a 
torque change with luminosity in this model. For 
4U 0142+61, as noted in Section |4~31 the magni- 
tude of the torque is decreasing while the pulsed 
flux is increasing. If propeller accretion is oc- 
curing, then the total flux should be decreasing; 
this can be checked (Gonzalez et al., in prepa- 
ration). We note that evidence against such a 
tor que /luminosity correla tion has been presented 
by (jGavriil fc Kaspill2004l ) for a different AXP. 

5. Summary 

Our continuing RXTE monitoring program has 
revealed a possibly new AXP variability phe- 
nomenon: 4U 0142+61 exhibited a slow but steady 
increase in its pulsed flux between 2002 May and 
2004 December, such that it has risen 36+3% over 
2.6 years in the 2—10 keV band. This is accom- 
panied by a softening of the spectrum (Gonza- 
lez et al. in preparation). Quasi-simultaneously, 
the pulse profile, which comprises two peaks hav- 
ing different spectra, has been evolving since 2000. 
In particular, the dip of emission between the two 
peaks has been rising since 2002, as if it is re- 
turning to its pre-2000 morphology in which there 
was no clear distinction between the peaks. The 
profile evolution translates to a reduction of the 
power in the Fourier harmonics of order higher 
than one since 2000. This is in contrast with 
the pulsed flux which seems to be moving away 
from the pre-2000 value. The evolution in the 
pulse profile is seen in the 2—4 keV band but not 
in the 6—8 keV band, presenting an interesting 
puzzle to the twisted magnetospher e mod el for 
magnetars. Intriguing orii et alJ I (1200,51 ) have 



suggested the pulsar suffered a glitch just before 
2000 on the basis of a single discrepant ASCA 
period measurement. Our phase-coherent timing 
using RXTE demonstrates that a glitch is plau- 
sible but not necessary to explain the data, but 
our pulse profile evolution analysis provides new 
evidence for such an event having occured. 

Physical interpretations that described well 
other observed long-term changes in AXP emis- 
sion (such as outburst afterglow or flux flares 
caused by an increased twist in the magneto- 
sphere) do not explain all the phenomena that 
we have observed. Most of our observations could 
be explained by the twisted magnetosphere model 
if the total flux of 4U 0142+61 is actually de- 
creasing. This would indicate a slow untwisting in 
the magnetosphere. Alternatively, if the total flux 
is increasing, a slow increase in the twist angle 
in the magnetosphere can account for the pulse 
profile simplification and for the post-Cycle 10 
events, but the changes in the spin-down rate and 
the softening of the spectrum would remain un- 
explained. Finally, the data could be explained 
by energy release following an energy deposition, 
perhaps due to a glitch that occured in the crust 
of the star sometime during the post Cycle 2 gap, 
although the energy deposited would have had to 
have been large. No matter what, the absence of 
profile evolution at high energies remains a puzzle. 
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Fig. 1. — Epochs of observations of 4U 0142+61 
with RXTE. 
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Year 
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Fig. 2. — Arrival time residuals for 4U 0142+61 
for the post-gap period, using the post-gap 
ephemeris given in Table [2 The residuals have 
RMS 2.3% of the pulse period. 
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Fig. 3. — Arrival time residuals for 4U 0142+61 
for all RXTE Cycles using the post-gap ephemeris. 
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Fig. 4. — Arrival time residuals for 4U 0142+61 
for all RXTE Cycles using the overall ephemeris 
(see Table [2]). The residuals have RMS 1.9% of 
the pulse period. 
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Fig. 5. — Pulse profiles in all RXTE Cycles in 
three different PCA energy bands (2—10 keV with 
64 bins across the profile, 2—4 keV with 64 bins 
across the profile, and 6—8 keV with 32 bins across 
the profile for clarity). Note that no observations 
were made in Cycle 4. Two full periods are shown. 
The normalization is such that the values of the 
lowest and highest bins in each profile are and 
1, respectively. The Cycle number is shown in the 
top right corner of each pulse profile plot. 
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Pulse Phase (cycles) Harmonic Number 

Fig. 6. — Fourier analysis of the pulse profiles 
in the 2—10 keV energy band. Left: Pulse pro- 
file curves made of the first six calculated Fourier 
components in each RXTE Cycle, superposed on 
the measured pulse profile points for that Cycle. 
Right: Harmonic content of the average pulse pro- 
files for each Cycle. Note that Cycle 3 was not 
included because of the poor signal-to-noise ratio. 
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Fig. 7. — Same as Figure [6] but for 2—4 keV. 
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Fig. 8. — Same as Figure [6] but for 6—8 keV. 
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Fig. 9. — Ratios of the Fourier amplitudes of the 
pulse profiles in two energy bands. Top: ratio 
of the Fourier amplitude of the second harmonic 
to that of the fundamental. Middle: ratio of the 
Fourier amplitude of the third harmonic to that 
of the fundamental. Bottom: ratio of the Fourier 
amplitude of the fourth harmonic to that of the 
fundamental. 
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Fig. 10. — Top: Pulsed flux evolution of 
4U 0142+61 in counts/s/PCU in the 2-10 kcV 
band. Each point corresponds to a full RXTE Cy- 
cle. Bottom: Pulsed flux evolution in erg/s/cm 2 
in the 2—10 keV band. See text for details. 
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Fig. 11. — Top: Superposed post-gap average 
pulse profiles in 2—10 keV (with six Fourier com- 
ponents included), scaled to give the appropriate 
average pulsed flux for each RXTE Cycle. Bot- 
tom: The evolution of the heights of three differ- 
ent features in the pulse profiles as a function of 
RXTE Cycle. The open triangles represent the 
maximum heights of the small peak in each Cycle. 
The open squares represent the maximum heights 
of the big peak. The filled circles represent the 
heights of the dip. 



25 




Fig. 12. — See caption for Figure [TT] but for 
2-4 keV. 
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Table 1 

Summary of RXTE Observations 



Observing 


Typical 


Typical 


Number of 


Total 


Cycle 


Exposure a 


Separation 21 


Obs. b 


Exposure First MJD — Last MJD d First Date — Last Date 




(ks) 


(weeks) 




(ks) 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 



11 
1 

20 

3 
7 
15 
5 
5 
5 



0.1 
4 



4 
14 
1 

15 
10 
9 
21 
27 
36 



45 
16 
20 


46 
54 
100 
124 
86 
120 



50170.469- 
50411.684- 
50893.083- 

51610.617- 
51986.347- 
52339.621- 
52726.197- 
53066.586- 
53438.151- 



50171.904 
50795.523 
50893.083 

51950.256 
52282.163 
52634.456 
53046.235 
53420.507 
53787.328 



29/03/1996-29/03/1996 
24/11/1996-13/12/1997 
21/03/1998-21/03/1998 



07/03/2000- 
18/03/2001- 
06/03/2002- 
28/03/2003- 
02/03/2004- 
09/03/2005- 



10/02/2001 
08/01/2002 
26/12/2002 
-11/02/2004 
19/02/2005 
21/02/2006 



a The exposure and separation are only approximate. 

b When the last digit of the observation ID of two successive data sets is different, the two data sets are considered separate 
observations. 

c The total exposure does not include Earth occultation periods. 

d First MJD and Last MJD are the epochs, in Modified Julian Days, of the first and the last observations in a Cycle. 
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Table 2 

Spin Parameters for 4U 0142+61 



Parameter 



Pre-Gap Ephemeris b Post-Gap Ephemeris 
Spanning Spanning 
Cycles 1 to 3 Cycles 5 to 10 



Possible Ephemeris 
Spanning 
All Cycles 



MJD range 
TO As 
v (Hz) 

v (1(T 14 Hz s- 1 ) 
v (1(T 23 Hz s" 2 ) 
d 3 v/dt 3 (10 -31 Hz s^ 3 ) 
d 4 v/dt 4 (10~ 39 Hz s- 4 ) 
d 5 u/dt 5 (10- 46 Hz s~ 5 ) 
d 6 u/dt 6 (10~ 54 Hz s~ 6 ) 
Epoch (MJD) 
RMS residual 



50170.693-50893.288 
19 

0.115099566(3) 
-2.659(3) 



50530.000000 
0.019 



51610.636-53787.372 
118 

0.1150969337(3) 
-2.6935(9) 
0.417(10) 



51704.000025 
0.023 



50170.693-53787.372 
137 

0.1150969304(2) 
-2.6514(7) 



51704.000000 
0.019 



-1.7(2) 
3.62(12) 

8.7(3) 
-5.01(13) 

6.6(4) 



a Numbers in parentheses are TEMPO-reported la uncertainties. 

b The pre-gap ephe meris reported here is slightly different from that reported in 
IGavriil fc KaspH (|2002h because here we take into account Cycle 1 and 3 observations. Note 
that both ephemerides return th e same number of pulsa r rotation cycles between the first and 
last pre-gap observations used by IGavriil fc Kaspil i|2002l ). 

c It is possible to find a different overall ephemeris after adding an arbitrary but constant time 
jump to all post-gap TOAs. 
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